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ABSTRACT 

Recent theoretical studies have shown that galaxies at high redshift are fed by cold, dense 
gas filaments, suggesting angular momentum transport by gas differs from that by dark mat- 
ter. Revisiting this issue using high-resolution cosmological hydrodynamics simulations with 
adaptive mesh refinement (AMR), we find that at the time of accretion, gas and dark matter do 
carry a similar amount of specific angular momentum, but that it is systematically higher than 
that of the dark matter halo as a whole. At high redshift, freshly accreted gas rapidly streams 
into the central region of the halo, directly depositing this large amount of angular momentum 
within a sphere of radius r = 0.1f? v ir- In contrast, dark matter particles pass through the cen- 
tral region unscathed, and a fraction of them ends up populating the outer regions of the halo 
( r/R v ir > 0.1), redistributing angular momentum in the process. As a result, large-scale mo- 
tions of the cosmic web have to be considered as the origin of gas angular momentum rather 
than its virialised dark matter halo host. This generic result holds for halos of all masses at all 
redshifts, as radiative cooling ensures that a significant fraction of baryons remain trapped at 
the centre of the halos. Despite this injection of angular momentum enriched gas, we predict 
an amount for stellar discs which is in fair agreement with observations at z=0. This arises 
because the total specific angular momentum of the baryons (gas and stars) remains close to 
that of dark matter halos. Indeed, our simulations indicate that any differential loss of angular 
momentum amplitude between the two components is minor even though dark matter halos 
continuously lose between half and two-thirds of their specific angular momentum modulus 
as they evolve. In light of our results, a substantial revision of the standard theory of disc 
formation seems to be required. We propose a new scenario where gas efficiently carries the 
angular momentum generated by large-scale structure motions deep inside dark matter halos, 
redistributing it only in the vicinity of the disc. 

Key words: galaxies: formation - galaxies: high-redshift - galaxies: spiral - galaxies: kine- 
matics and dynamics - galaxies: fundamental parameters - cosmology: theory 


1 INTRODUCTION 

Along with mass, angular momentum growth shapes the funda- 
mental properties of galaxies. Using linear theory, one can show 
that density perturbations acquire angular momentum through their 
asymmetric interactions with the large r scale tidal field jPeeblesI 
1 19691 : lDoroshkevicf]| 197d : [whitdl 1 9841) . Thus one expects the an- 
gular momentum of the Lagrangian region encompassing a fu- 
ture virialized structure to scale like a 2 (t)D(t) until it decou- 
ples from the Hubble expansion, where a and D are the scale 
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factor of the universe and the linear growth rate of density per- 
turbations respectively (IWhitel 179841 ). Further assuming that the 
gas experiences the same tidal field as the dark matter, they 
should initially share their specific angular momentum. Provided 
this specific angular momentum is conserved as the gas radia- 
tively cools, early studies were able to reasonably match crucial 
observed relations, such as the Tully-Fisher relation or the size- 
rotation velocity relat ion dFall & Efstathioulll980l : lDalcanton et al.l 
1 19971: Mo et al.l 1 1 9981) . In spite of sev eral sever e short comings in 
the theory pointed out by authors like iHoffmanl dl986l) (including 
no account of secondary infall and mergers between virialized ob- 
jects), this success encouraged all (semi-analytic) galaxy forma- 
tion models published to date (e.g. ICole et aLll200d : lHatton et al.l 
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2003; |Croton_etai]_ 2006^ Monacoetal 20071; Somerville et al .1 


2008 ; iDutton & van den Boschl 20091 ; Khochfar & Silk! 201 lh to 
rely on the core assumption that gas and dark matter contained 
within the same virialized structure split specific angular momen- 
tum equally. However, the recent findings that the vast majority of 
galaxies are mainly fed gas by cold, thin and dense filaments which 
penetrate deep inside the virial radius of their host dark matter 
halo fe.g. IKeres et al.|[20t)I bcvirk et al.ll2008l: iDekel et al.ll2009l: 


IPowell et al.ll201 lh begs for a reassessment of the validity of the 
classical angular momentum scenario. 

According to the standard picture of galaxy formation, gas 
is first shock-heated to the virial temperature of its dark matter 
halo host as this latter collapses, and subsequent accreted mate- 
rial encounters the resulting viria l shock as it penetrat es the halo 
dRees & Ostrikerll 1 9771: fsilMl 1 9771: 1 White & Reesll 1 9781) . Since the 
shock-heated gas is more or less spatially re-distributed as the dark 
matter before it can cool, one is led to logically postulate that its 
specific angular momentum closely tracks that of the dark matter 
halo. Non-radiative hydrodynamic s cosmological simulations have 
confirmed that this is the case fe.g. Ivan den Bosch et al.ll2002l) . 

However, as soon as the gas is permitted to radiatively cool 
and form stars, this tight correlation between gas and dark matter 
spin should break down. For instance, a fraction of gas accreted at 
earlier times with lower angular momentum will be converted into 
stars, elevating the specific angular momentum of the remaining 
gas button & van den Boschl 12009 ^). Meanwhile the specific an- 
gular momentum of the dark matter halo is expected to either be 
frozen in after t urn around or decline due to the p resence of exter- 
nal torques (c.f. |Peebles|[~1969l : iBook et aT]|20 1 II) . Hence it is not 
very surprising that radiative hydrodynamics cosmological simula- 
tions show a difference in the spin p arameters between gas and dark 
matter. Indeed. [Stewart et~ahl 1 20 1 lh have recently reported that for 
a pair of resimulated halos, cold gas had 3 to 5jime_s more specific 
angular momentum than dark matter (see also lSales et al . 201Gh . In 
this context, the questions which naturally arise are: (i) what are the 
mechanisms that segregate dark matter and gas angular momenta? 
(ii) is this segregation universal or does it depend on halo proper- 
ties? and (iii) does it evolve with redshift? Since the gas is the main 
agent for transferring angular momentum from (super) halo scales 
down to central galaxies, it is key to investigate the evolution of its 
angular momentum in the region extending between the galaxy and 
the virial radius of the dark matter halo (i.e. 0.17? v ir < r < R v i r ). 

Within galaxy size virialised halos, the picture of gas ac- 
cretion has recently been significantly revi sed, with the redis- 
covery of an idea first put forward by iBinnevI (1 1977h that in- 
falling gas is never shock-heated to the virial temperature but in- 
stead flow s through an ‘isothermal’ shock, reaching the galactic 
disc cold dKat z et al.ll2003l: iBimboim & Dekell 12003 ; IKeres et al.l 
l2005l : l0cvirk et all2008l : Brooks et al.ll2009lk Moreover, this gas is 
brought along anisotropic narrow streams which persist deep within 
the halo, contrary to their broader dark ma tter cou n terpar ts which 
rapidly fade away at the virial radius ll Powell et al.ll201 fl) . There- 
fore, one wonders how differently this filamentary gas advects an- 
gular momentum throughout the halo. 

With such questions in mind, we use a set of high-resolution 
cosmological simulations to revisit two key assumptions of disc 
formation theory, namely that gas carries the same amount of spe- 
cific angular momentum as its host dark matter halo and that this 
specific angular momentum is conserved as this material gets ac- 
creted onto the central disc. The details of the simulations are de- 
scribed in Section 2. In Section 3, we present the time evolution 
of angular momentum of baryons and non-baryonic matter along 



Figure 1. Projected densities of the dark matter (top) and gas (bottom) cen- 
tred on a halo of ~ 1 0 1 1 A7q at z = 3 from the NUT simulation with super- 
nova feedback (NutFB). A 6R v ir region is projected and a figure with as- 
pect ratio 1:2 is displayed, where the longer length corresponds to = 6i? v i r 
(= 224 kpc). At z = 3 the scale-height of the disc as well as the large scale 
filamentary structure is well resolved in the simulation. It can be seen that 
dark matter filaments are broader than gas filaments. The gas filaments are 
not destroyed by supernova explosions. 


with a comparison to observational data at 2 = 0. We also dis- 
cuss in this section how our results depend on halo properties and 
evolve with redshift. Finally, we discuss the discrepancy between 
dark matter and gas specific angular momenta in Section 4 and con- 
clude in Section 5. 


2 SIMULATIONS 


We use the Eulerian AMR code, Ramses (lTevssieill2002l) . to inves- 
tigate the angular momentum evolution of baryons and dark matter. 
Ramses uses a second-order Godunov scheme to solve the Euler 
equations, and an adaptive particle-mesh method to solve the Pois- 
son equation. Since the outcome of hydrodynamics simulations is 
subject to resolution and physical ingredients, we make use of five 
simulations to draw robu st conclusio ns. These simulations include 
the Nut se ries i PowelletaL 201 lh. the HoRizoN-Ma reNostrum 
simulation llOcvirk et al. boQS : Devriendt et al.| l2010h . and two 
other cosmological simulations (Cosmo25, Cosmo50, Dubois et 
al. in prep.) of fairly large volumes but with lower resolution than 
the HoRizoN-MareNostrum simulation. Whereas simulations in the 
Nut series focus on the evolution of an individual. Milky Way-like 
galaxy, the HoRizoN-MareNostrum, Cosmo25 and Cosmo50 sim- 
ulations give the statistical properties of galaxies spanning a wide 
range of halo mass. Physical ingredients common to most of our 
simulations include reionisation, cooling, star formation and super- 
nova feedback. These are summarised in Table 1, along with simu- 
lation parameters. 

To model the reionisation of the Universe, a uniform UV 
background radiation field (modelled as a heating term in the en- 
ergy equa tion) is turned on at high redshift (see Table [TJ follow- 
ing lHaardt & Madaul (1 1 9961). Gas dissipates energy t hrough atomic 
cooling down to 10 4 K i Sutherland & Doniu]l 1 9931) . For the Nut 
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simulations metal line cooling can lower the gas temperature fur- 
ther (below ~ 10 K). When the gas density in a grid cell exceeds 
a given threshold (nth, see Table 1), star particles are spawned 
by a Poisson process according to a Schmidt law with a 1 % ef- 
ficiency of the star form ation per free-fall time jKennicutdl 1 9981 : 
iDubois & Tevssie3l2008l) . This threshold density is chosen so that 
it is inferior or equal to the maximal Jeans density reached on the 
finest level. In the simulations with supernova feedback, after ~ 10 
Myrs, massive stars undergo Type II supernova explosions, releas- 
ing half of their 10 51 ergs into their surr oundings as kinetic energ y 
and the other half as thermal energy iDubois & Tevssieg|2008l ). 
During this phase, processed heavy elements are dispersed, enrich- 
ing the interstellar and intergalactic medium. In what follows, we 
elaborate on the details of each simulation. 

The Nut series make use of the resimulation (also called 
‘zoom’) technique to follow the evo luti on of a M ilky Way-type 
galaxy in a ACDM cosmology. IPowell et alj 120 1 ll) reported on 
high redshift results (z > 9) from the ultra-high resolution resimula- 
tions in the Nut series. These reached a maximum physical spatial 
resolution of ~ 0.5 pc at all times. To track the evolution of the 
galaxy down to lower redshifts, the Nut suite also includes three 
resimulations with lower spatial resolution (but identical DM par- 
ticle mass resolution) and the following physics: (i) adiabatic with a 
uniform UV background turned on instantaneously at z= 8.5 (Nu- 
tAD) (ii) cooling, star formation, and UV background (NutCO) and 
(iii) same as (ii) but with supernova feedback and metal enrichment 
(NutFB). NutAD and NutCO have maximum 48 pc (physical) res- 
olution at all times and reach z = 0, whereas NutFB has maximum 
12 pc (physical) resolution and only reaches z = 3. Fig.[T]shows a 
snapshot of the DM (upper panel) and gas (lower panel) of a region 
centred on a halo hosting a Milky Way-type galaxy at z = 3 in 
NutFB. 

We recall the important details of the Nut resimulations 
here. The simulation volume is a 9/i -1 Mpc com oving pe ri odic 
box e volving according to a WMAP5 cosmology JPunklev et al.1 
1 20091) (a s = 0.8, D m = 0.258, Da = 0.742, h = 

.ffo/(100kms -1 Mpc -1 J = 0.72). Initi al conditions are gener- 
ated using MPgraficj Pranetet al]2008l) . a parallel version of the 
Grafic package jBertschingej|20oH )7 Within this volume we iden- 
tify the region where a Milky Way-like galaxy (halo virial mass 
of ~ 5 x 10 11 Mq at a = 0) will form. This region encompasses 
a volume of side length ~ 2.7 IP 1 Mpc. While the root grid for 
the entire simulation volume is 128 3 , within the (~2.7 h -1 Mpc) 3 
region, we place an additional three nested grids, giving an equiva- 
lent resolution of 1024 3 dark matter particles, each with mass Mdm 
~ 5 x 10 4 Mq. To fix the maximum physical resolution to a con- 
stant value (12 pc for NutFB and 48 pc for NutCO and NutAD) as 
the universe expands and the simulation evolves, we further refine 
on the finest fixed grid within the 2.7 IP 1 Mpc 3 region according 
to a quasi-Lagrangian strategy, i.e. when the number of dark matter 
particles in a cell reaches 8 or equivalently when the baryon plus 
dark matter density in a cell increases by a factor of 8. Table[I]lists 
the maximum level triggered for each simulation. Because of the 
higher spatial resolution in NutFB it uses a higher density thresh- 
old for star formation (nn,th = 400 cm -3 ) than the NutCO run 
(riH.th = 10 cm -3 ). The star particle mass (~ 2 — 3 x 10 4 Mq) is 
determined by the combinati on of minimum grid size and density 
threshold for star formation ( IDubois & Tevssienl200a) . In NutFB, 
we assume that every supernova bubble with an initial radius of 32 
pc sweeps up the same amount of gas as that initially locked in the 
star particles. This is usually expressed as a mass loading factor of 
unity ( rj = 1). 


The other three simulations are large volume cosmological 
simulations, and as such are performed with lower spatial (1-2 
kpc) and mass resolution. More specifically, the mass of each dark 
matter particle is jtidm — 10 7 Mq for the HoRizoN-MareNostrum 
simulation, 8 x 10 7 Mq for the intermediate size run (Cosmo25) 
and 6 x 10 8 Mq for the Cosmo50 run. As in the NutFB run. su- 
pernova feedback and UV background heating are included in the 
simulations, but the radius of the initial supernova bubble is set to 
twice the minimum size of the grid (see Table |T}. Note also that 
the adopted cosmology for the HoRizoN-MareNostrum simulation 
(WMAP1) is different from the others, but as we will show this has 
very little impact on our resul ts, if at all. We refer in terested read- 
ers to lOcvirk et al.l J2008t) and lDevriendt et alj J 20 1 (ll) for a detailed 
description of the HoRizoN-MareNostrum simulation set-up. 

In all the simula tions, we identify (sub) haloes using the Adap- 
taHop algorithm lAubert et al.ll2004l ). which is based on the detec- 
tion of peaks and saddle points in the dark matter density field, 
sup plemented b y the most-massive subhalo algorithm developed 
by [Tweed et alj ( I2009f) . The virial radius of halos is defined as the 
maximal radius within which the virial theorem is satisfied to bet- 
ter than 20%. We further define gas belonging to a satellite galaxy 
as gas residing within the half-mass radius of its host DM satel- 
lite halo. The centre of a halo, which we use to compute angular 
momentum, is defined as the centre of mass of dark matter and 
baryons. The mean motion of the halos is determined by comput- 
ing the centre of mass velocity of dark matter particles, gas, and 
stars within their virial radii. 


3 THE ANGULAR MOMENTUM OF VIRIALIZED 
HALOS 

To understand how galaxies acquire their angular momentum, we 
study the angular momentum evolution of the different components 
(gas, dark matter, stars) inside their host halos. We compute the 
specific angular momentum vectors as: 

r Ei m * f -xvi 

Ei m i 

where r,: is the radial distance from the centre of mass of the halo 
(includes dark matter and baryons), Vj is the peculiar velocity and 
Trii is the mass of the i-th dark matter (star) particle or gas cell. 
In what follows we use different subscripts to denote the specific 
angular momentum of different components in different regions of 
the halo. These are summarised in Table 2. 


3.1 Adiabaticity and the cosmic origin of angular momentum 

We begin our investigation of the evolution of angular momentum 
of the various components of a virialized halo with the simplest 
(in terms of physics) high resolution NutAD run. Since gas cannot 
cool radiatively in this simulation, newly accreted material is shock 
heated by the pressure supported intra-halo medium. As a result, 
its radially oriented initial velocity is isotropized, and drives the 
gas density field towards spherical symmetry. Accreted dark matter 
particles are also more or less isotropically redistributed within the 
host halo by the collisionless violent relaxation process. Given that 
(i) gas and dark matter within the halo experience the same larger- 
scale torques (e.g. [Peebles! 1969i: lBook et alj201ll) . and (ii) gas and 
dark matter are driven by the gravitational collapse to a very sim- 
ilar equilibrium distribution (i.e. to a good approximation that of 
an isothermal sphere since the total amount of angular momentum 
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Table 1. Summary of simulation parameters and physical ingredients. From left to right, columns are as follows: simulation name, comoving box size, number 
of dark matter particles, level of the root grid, level of the finest grid, minimum grid size, dark matter particle mass, star particle mass, threshold density for 
star formation, redshift down to which the simulation is carried out, indication of whether supernova feedback is included, redshift at which UV background 
heating is initiated, and remarks. 


Simulations 

L 

(Mpc/h) 

Adm 

^min 

^max 

Aa-min 

(pc) 

m DM 

(M&) 

f'ft-star 

(M q ) 

™th 

(H/cm 3 ) 

^end 

SN 

UV 

Remarks 

NutFB 

9 

1024 3 

7 

20 

12 

5 x 10 4 

2 x 10 4 

400 

3 

Y 

N 

II 

OO 

Ui 

Resimulation 

NutCO 

9 

1024 3 

7 

18 

48 

5 x 10 4 

3 x 10 4 

10 

0 

- 

z=8.5 

Resimulation 

NutAD 

9 

1024 3 

7 

18 

48 

5 x 10 4 

- 

- 

0 

- 

z=8.5 

Resimulation 

MareNostrum 

50 

1024 3 

10 

16 

1090 

1 x 10 7 

2 x 10 6 

0.1 

1.5 

Y 

z=8.5 


Cosmo25 

25 

256 3 

8 

15 

1090 

8 x 10 7 

4 x 10 6 

0.1 

0 

Y 

z=10.5 


Cosmo50 

50 

256 3 

8 

15 

2180 

6 x 10 s 

3 x 10 7 

0.1 

0 

Y 

z=10.5 



Table 2. Notations (first column) for specific angular momenta used in this 
study. The first subscript indicates the component used for the calculation 
(DM, gas, stars, baryons). The second column specifies the region (nor- 
malised to the halo virial radius) over which the specific angular momentum 
is calculated. The default region over which we measure specific angular 
momentum is a sphere extending from the centre of mass of the halo to the 
virial radius. A second subscript indicates whether we exclude 0 gas, out) or 
include (jgas,cen) the gas in the central region {r / R w i r < 0.1) of the halo 
for the specific angular momentum calculation. The third column lists the 
components included in the specific angular momentum calculation. 


Notation 

Spatial extent 

Component 

2 dm 

r/Rv ir < 1 

dark matter 

J gas 

r/Ror < 1 

gas 

2gas,out 

0.1 < r/H vir < 1 

gas 

2gas,cen 

r / R v i r < 0.1 

gas 

J star 

r/Rv ir < 1 

all stars except those in satellite galaxies 

2bar 

r/Rv ir < 1 

all stars (satellite galaxies included) + gas 


provided by tidal torques is very limited), we expect them to have 
similar j. Indeed, Fig. [2] shows that jgas,out closely tracks jdm re- 
gardless of whether or not mergers, easily identified by large jumps 
in j, occur. When the dark matter angular momentum is expressed 
in terms of the spin parameter (A' = rim /v^^irT^. lBullock et al] 
hooih , we recover the typical value of ~ 0.04 at all times with 
fluctuations of up to a factor of ~ 2 around this value. 

An interesting feature in Fig.[2]is that jd m increases with time 
(albeit at a reduced rate below z = 3), implying that late infall car- 
ries a larger amount of angular momentum. This is not a completely 
unexpected result given that dark matter haloes are kn own to expe- 
rience little evolution of spin parameter with time (e.g. lPeirani et al .1 
|2004 As their mass and radius grow, so must their angular mo- 
mentum. However, it is not trivial to understand why the late accre- 
tion has larger j. A naive answer is that material with larger angu- 
lar momentum takes more time to reach the potential well, but this 
does not explain how it acquired such a large angular momentum 
in the first place. We discuss in detail the cosmic evolution of the 
angular momentum of halos (which ne cessarily goe s beyond tidal 
torque theory) in a companion paper IPichon et al.ll20 1 1J), but for 
completeness’ sake, we briefly outline the main idea in the para- 
graph below. 

The dynamics of the gas and dark matter flowing along what 
has been dubbed the ‘cosmic web’ can be understood as the 
anisotropic time evolution of the initial gaussian random gravi- 
tational field. Cosmological struct ures originate as peaks in the 
associated initial density field ( IBond et al.ll 1 9960 . which are con- 


look-back time [Gyrl 

4 8 10 12 13 



Figure 2. Evolution of the specific angular momentum (j) of a dark matter 
halo hosting a Milky Way-type galaxy in the adiabatic run (NutAD). Dif- 
ferent colours denote different components: dark matter (y'dm. thick grey), 
halo gas O' ga s,out, 0.1 < r/f? v ; r < 1, thick blue). Also plotted as a blue 
dashed line is j of all the gas inside the virial radius, y'gas. We also dis- 
play the specific angular momentum corresponding to a reduced spin value 
A' = 0.04 (red dashed line) for the dark matter halo. The specific angular 
momentum of the halo gas O'gas.out) closely follows that of the dark matter 
halo Odm) in ibis non-radiative simulation. 


nected to other peaks through peak patches. The bulk motion of a 
peak patch is determined by the gradient of the large-scale poten- 
tial. The latter also drives the motions of filaments which exist at 
the intersection of at least three void patches and connect peaks. 
As a consequence of the asymmetry between voids, the gas and 
dark matter flowing out of these voids acquire a transverse veloc- 
ity when they intersect at a filament. This transverse velocity is the 
seed of a halo’s angular momentum which is then advected along 
the filaments all the way into the halo sitting on the peak. Since 
the transverse velocity along a filament is constant to zeroth or- 
der approximation, the material initially located further away from 
the peak will naturally contribute more angular momentum. Note 
that the infall of matter along such filaments will coherently con- 
tribute an increasing amount of angular momentum as time goes on 
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Figure 3. Left panel: evolution of the specific angular momentum modulus ( j ) from the NutFB run. Different colours denote different components: dark matter 
Odm’ thick grey), gas at 0.1 < r/R v i r < 1 Ogas.out, thick blue), central gas (r/R v i r < 0.1, j gas ,cen, thin blue), stars (jst ar , thin orange), all baryons 
in the halo Obar’ dot-dash red). Note that the specific angular momentum calculations for dark matter Odm)’ gas Ogas.out) and baryons Obar) includes the 
contribution from satellite galaxies, while it is excluded for the measurement of j B t ar and Jgas.cen. Angular momentum corresponding to a constant reduced 
spin value of A' = 0.04 is included as a grey dashed line. The specific angular momentum of gas Ogas.out) is larger than that of its host dark matter halo 
Odm). suggesting that the angular momentum of the gas is not acquired during the collapse of this latter. Right panel: specific angular momentum along the 
spin axis of the gas in the central region Ogas.cen). By definition, j ■ y'gas.cen for gas in the central region in this panel (thin blue line) is the same as j gas ,cen 
in the left panel. We denote negative values as dashed lines. Different angular momentum amplitudes between lines of the same type in the left and right panels 
indicate the misalignment between the central gas and each component. 


(shown in Fig. [2}, since the filament preserves its orientation over 
long stretches of time. 


3.2 Radiative cooling, supernovae feedback, and the 
evolution of angular momentum at high redshift 

In the real Universe, gas can radiatively cool and form stars, which 
should therefore be considered as a sink for gas mass and angular 
momen tum. Moreove r, at high redshift, star formation is very active 
dBouwe ns et al .120091) and the effect of supernova explosions on the 
evolution of the halo gas may prove important. Thus, we make use 
of the high-resolution NutFB run which includes these processes 
to investigate the evolution of angular momentum of a DM halo at 
high redshift ( z > 3). 

The left panel of Fig. [3] shows that j ga s ,out and jdm increase 
steadily, although with occasional rapid variations located around 
look-back times (fib) of fib =12.9, 12.4, and 11.7 Gyr. These an- 
gular momentum ‘jumps’ are generated by episodes of enhanced 
accr etion of dense gas contained in satellite galaxies and filaments 
(e.g. lMaller et alJl2002l : |Peirani et alj|2004l) . Indeed, the most sig- 
nificant halo merger occurs at fib — 11.7 Gyr (satellite to host 
halo mass ratio Ms a t/Mhost ~ 0.1) and the majority of the other 
12 minor mergers (0.01 < M sa t/Mhost < 0.1) take place be- 
tween 12.5 )$ fib ^ 13.0 Gyr. In addition to these mergers, mate- 
rial accreted along filaments also contributes to the rapid variation 
of 7 since the density in the gas filaments varies (see also Fig. 2 of 
IBrooks et alj20Qp) . 

However, the key feature of Fig.[3](left panel) is that y ga s ,out, 
although tightly correlated with jdm, is always larger by a factor 


of a few (see also QtenetaL 20031 Sharaia^^temmetz 2005|; 


iGottlober & Yepesl 120071 : ISales et alj 20 id : IStewart et all 201 lh . 
This is in contradiction with the picture in which the angular mo- 
mentum of the halo gas is acquired during the virialisation of its 
host dark matter halo and suggests that if the gas and the DM halo 
interact at all and exchange angular momentum during the collapse, 
the effect is likely to be unimportant, at least in the outer regions 
of the halo. Indeed this is more in line with the picture according 
to which gas accreted at the virial radius supe rsonically streams 
(almost) directly i nto the central galaxy (e.g. IBrooks et al .11 20091 : 
tPowell et al JI201 lh . Moreover, Fig. [3] (left panel) also shows that 
jg as, ceii (thin blue solid line) is systematically smaller than y ga s.out 
(thick blue solid line). There are two possible explanations for this. 
First, the specific angular momentum modulus of the gas is not con- 
served during its gravitational collapse and subsequent accretion 
onto the c entral galaxy because .j gas , ou t is misaligned w i th.j gas .cen 
(also see Ivan den Bosch et aD 120021: lHahn et al] I20T0I: iBett et al.l 
bold) . The second possibility is that unlike the gas in the outer part 
of the halo which yields an instantaneous picture of the accreted 
angular momentum, the central gas reflects the angular momentum 
accretion history of the gas. We now discuss each of these possibil- 
ities in turn. 


The difference in angular momentum alignment between outer 
and inner gas is indicated in the right panel of Fig. [3] where we 
plot the component of j gas ,out (thick solid blue line) parallel to 
j g as,cen • This figure clearly shows that this component is remark- 
ably well preserved (as one would expect for an isolated axisym- 
metric system) and it is therefore tempting to conclude that mis- 
alignment accounts for most of the difference between j gas , ou t 
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and jgas.cen- However, the time evolution of the projected com- 
ponent of j gas, out along jgas,cen direction does not correlate very 
well with the time evolution of j g as ,cen. This behaviour is natu- 
rally interpreted as a gas accretion history effect. Indeed, as the 
mass of gas in the low amplitude specific angular momentum cen- 
tral region (r/R v i r < 0.1) is comparable (between a factor 2-3 
smaller) to that in the outer reservoir (0.1 < r/R v i r < 1) (see 
Appendix for details on the evolution of the gas mass in these re- 
gions), even if jgas.out was perfectly aligned with j ga s,cen, the re- 
sulting modulus of the specific angular momentum of the inner re- 
gion after accretion would naturally be reduced compared to that 
of the freshly accreted gas. As a result, notwithstanding mergers, 
igaa.cen grows more smoothly than j s as , ou t and takes more time 
to reach a given level of specific angular momentum. This inte- 
grated effect is exacerbated for the stellar component, which dis- 
plays an even smaller amount of specific angular momentum than 
Jgas.cen- Since the global star formation timescale is much longer 
than t he dynamical timescale of the central region (e.g.. lKennicuttl 
1 19981) . most of the stars will form from older gas, with a lower 
angular momentum amplitude. Even though newly formed stars 
carry more angular momentum, they represent a smaller fraction 
of the total stellar mass than the central gas mass compared to 
that of the total halo gas. Therefore, a supplementary time delay 
to reach a given level of j and a smoother growth in specific an- 
g ular moment u m amplitude of the stellar component is inevitable 
button & van den Boschj|2009l) . 

Finally, when we measure the specific angular momentum of 
the total baryonic component, i.e. gas and stars, ( j'bar : thin red dot- 
dashed line in the left panel of Fig. [3}, we find that its amplitude 
most closely tracks that of the specific angular momentum of the 
DM (jdm, thick grey solid line) which we know reflects the mass 
assembly history of the halo. We are not implying that |jdm| of 
individual dark matter particles is conserved during the collapse, 
but simply that, as the baryon angular momentum is advected deep 
within the halo, we can only probe the accretion history of the 
gas by including all baryons within the inner region. When this 
effect is accounted for, there is no indication in our simulation 
that baryons lose more (or less) specific angular momentum am- 
plitude than DM during gravitational collapse of the halo. Interest- 
ingly, the specific angular momentum of the central gas has a value 
close to that of dark matter for virialise d halos typically measured 
in N-body simulations (A' = 0.04) (e.g. lEfstathiou & Jonejl 1 97 c j ; 
IBames & Efstathioulll987r) . It should be noted, however, that it is 
established with the aid of gas within 0.1 < r/_R v ir < 1, which 
has a larger j than the dark matter. We therefore surmise that if this 
gas had started with the same specific angular momentum as the 
dark matter halo to begin with, the central gas would have been left 
with roughly half of the typical DM specific angular momentum. 
We come back to the crucial issue of specific angular momentum 
segregation between gas and DM in Section[4]but now address the 
question of the importance of feedback in regulating j'gas.cen- 

Using a couple of cosmological hydrodynamics resimulations, 
one of a bulge-domin ated galaxy a n d ano ther of a disk-dominated 
galaxy at redshift 0, IZavala et al J 120081 ) show that the specific 
angular momentum of the disc galaxy closely follows that of its 
dark matter halo. They argue that this is because the stellar feed- 
back in the simulation which produces the disk-dominated galaxy 
is strong enough to heat the gas within the halo thereby prevent- 
ing it from cooling too much at early times, fragmenting into sub- 
galactic clumps, and transferring its angular momentum to the dark 
matter halo by dynamical friction. Instead this hot gas accumulates 
in the halo where it acquires the specific angular momentum of 



look-back time [Gyr] 

Figure 4. Fraction of total angular momentum of gas at 0.1 < r / R vlr < 1 
from various components. We divide the gas into four different phases: i) 
cold (T < 10 5 K), dense filamentary (n H fy < njj < n H t j,), ii) cold 
diffuse (rtn < rr k . fi [ ), hi) hot (T > 10 5 K) dense, and iv) hot diffuse 
component (see Equation m for the definition of njj m). Note that these 
quantities do not include the gas belonging to satellite galaxies. Grey solid 
line corresponds to the total fraction of angular momentum of gas belonging 
to the host halo. Cold, dense filamentary accretion accounts for most of the 
angular momentum in the halo. 


the dark matter before it eventually cools and re-collapses onto the 
central galactic disc, all the while conserving its angular momen- 
tum. However, in the first cosmological resimulation where indi- 
vidual supernovae remnants are resolve d, very littlegas is heated 
by stellar feedback at high redshift (I Powell et al.ll2011l) . Moreover, 
these authors show th at the cold filaments which supply most of 
the gas to the galaxy llKeres et al.ll20 05ll2009l : iDekel & Bimboiml 
120061: lOcvirk et al.ll2008l: iBrooks et al.ll2009l) are not disrupted by 
su pernova feedback. This r esult h as been independently confirmed 
by iFaucher- Giguere et alj d20 1ll) who use d the constant velocity 
wind model of Springel & Hemquistl d2003l) to model feedback and 
found that at high redshift ( z > 3) their supernova-driven galac- 
tic winds with low mass loading (17 = 1) were not able to sup- 
press the accretion of cold, dense gas in intermediate-mass haloes 
(M v i r ~ 10 11 - 12 Mg). Only when they assumed extreme param- 
eters for their wind (i.e. loaded them with double the mass that is 
turned into stars, and drove them with double the amount of en- 
ergy than that available from their supernova explosions) were they 
able to shut down cold, dense accretion. However this took place 
at the expense of the baryonic mass function which_ended up dra- 
matically undershooting the observations of e.g. iBell et ali 12003 ) 
at z = 0. Purely based on energetic grounds, feedback from ac- 
tive galactic nu clei m ay be able to blast the filamentary structure 
ivan de Voort et alj20ld) . but it is unclear how geometrical effects 
will affect their capacity to do so. A collimated jet for instance 
will not, in general, deposit enough of this energy in the vicinity of 
the filaments because of their small covering factor (a few percent, 
iKimm et al JI201 ll) . In any case, i n the simulat i on pr esented here, 
and contrary to the conclusions of lZavala et al .1 d2008h . the specific 
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angular momentum of the gas in the galaxy (thin blue line in the 
left panel of Fig. El is similar (slightly larger) to that in its DM halo 
host, even though supernova feedback is moderate. This raises the 
question of the importance of the numerical technique employed 
to assess the fragmentation and transfer of angular momentum in 
disks. Although this is beyond the scope of this paper, as we are 
chiefly interested in the outer parts of the halo, we believe that our 
simulation provides a more correct answer simply because we have 
better resolution. We ref e r the interested reader to the recent pa- 
per bv lCommercon et~ahl J2008ll for a thorough comparative study 
between SPH and AMR as to how a high artificial viscosity and 
numerical noise undermine the capacity of the SPH technique to 
conserve angular momentum of disks when the number of particles 
(total and in the kernel) is not high enough. 

Since accretion through cold , dense filamentary gas is thought 
to do minate at high redshift (e.g. iBrooks et al JI2009I ; I Powell et al.l 
1201 lh . the larger / gas ,out relative to /dm is also expected to be as- 
sociated with the dense filamentary component. To test this idea, as 
well as assess the importance of supernova feedback, we divide the 
gas into four different phases according to temperature and density 
as 

(i) cold dense ( T < 10 5 K and riH.fli < riH < nth) 

(ii) cold diffuse (T < 10 5 Kand hh < n.H,fu) 

(iii) hot dense (T > 10 S K and nn.fli < riH < n t h) 

(iv) hot diffuse (T > 10 S K and hh < riH.fii). 

where we define the cold dense filamentary structure using the 
nn.fii < riH < riH.th gas density cut. The lower density bound 
(nH,m) is chosen as 

nH.fli = SipfhaiXn/mn, (2) 

where /bar, Xu, and mH are the universal baryon fraction 
(flb/fim), the primordial hydrogen mass fraction and the mass 
of the hydrogen atom, respectively. The dependence of the lower 
bound of the filament density on the background density (p(z)) 
is motivated by the fact that the filamentary structure acquires its 
properties on large scales which are not gravitationally decoupled 
from the expansion of the Universe. The parameter 5f determines 
the overdensity of the filamentary gas. We find that <Sf ~ 100 
reasonably identifies the filamentary structure in the Nut simula- 
tions, which corresponds to hydrogen number densities higher than 
riH — 0.02 and ~ 0.001 cm~ 3 at z=9 and z=3, respectively. Note 
that, although this is not a worry in the Nut simulations because 
we spatially resolve the filaments, 5f will generally be resolution- 
dependent in cosmological simulations. In a rather obvious fashion, 
a lack of spatial resolution will artificially broaden the filamentary 
structure and thus reduce its density as mass needs to be conserved. 

Fig.[4]shows the contributions from the four different phases to 
the total angular momentum of the gas located at 0.1 < r/R v i r < 
1. We emphasize that gas belonging to satellite galaxies is excluded 
from the measurements. It is clear from this figure that cold fila- 
mentary accretion is primarily responsible for the larger / gas , out. 
The cold filamentary gas alone accounts for more than half of the 
gas angular momentum in the region. Thus, it can safely be con- 
cluded that the cold, dense filamentary accretion carries material 
with larger j than the dark matter halo from outside the virial ra- 
dius to the inner region of haloes without the gas having much time 
at all to interact with the dark matter. Gas belonging to satellite 
galaxies occasionally contributes a significant fraction of J (dif- 
ference between the solid grey line and a horizontal line drawn 
for J • /gas, out / J g as,out = 1 in Fig [4]). However, its impact on 
the specific angular momentum is less significant than that of the 


cold, dense filamentary accretion at z > 3 for two reasons. First, 
whilst cold, dense filamentary gas flows into the central gaseous 
disc rapidly, gas gravitationally bound to a satellite galaxy orbits 
around it for the time it takes dynamical friction to drag the satellite 
galaxy down. As a consequence of this process, satellite gas angular 
momentum in Fig.[4]cannot be directly converted into the actual an- 
gular momentum of gas accreted onto the central galaxy. Second, 
satellites are accreted along the gas filaments, and hence part of 
the gas which we conservatively assigned to these satellite galaxies 
could be regarded as filamentary gas. In this sense, the estimate of J 
for the satellites in Fig.[4]should be considered as an upper bound, 
and thus the actual fraction of the total angular momentum con- 
tributed by cold filamentary accretion may be slightly higher than 
reported in this work. Finally, the contribution to the total angular 
momentum from other phases (hot diffuse, hot dense, and cold dif- 
fuse) to J turns out to be minor for a DM halo of this mass in this 
redshift range. Whilst the lack of angular momentum of the cold 
diffuse phase originates from its almost isotropic accretion which 
cancels its angular momentum to a high degree, that of the hot gas 
phase is mainly caused by the small mass fraction of hot gas. This 
result is not very surprising because our DM hal o, at 1 0 11 is not 
massive enough to susta in a stable virial shock JBirnboim & Dekell 
120031 : iKeres et aPl2005h . and supernova explosions in simulations 
are notoriou sly inefficient at ejecting l arge amounts of hot gas 
in the halo 1 Mac Low & Ferraral 1999^ Dubois & Tevssieil 120081 : 


IPowell et all201 it ' Faucher-Giguere et al.ll201 1 ). 


3.3 Shock heating and the evolution of angular momentum 
at low redshift 

Gas infall at low r edshift is dominated b y smooth accretion rather 
than mergers (e. g. iFakh ouri et al]|201Cl and references therein). 
As discussed by IPichon et al. r honi) . this substantially increases 
the specific angular momentum of gas in the halo by delivering 
it coherently via the large-scale cosmic web. This is illustrated in 
Fig. El where the time evolution of / gas , 0 ut for the Milky Way- 
like halo in the NutCcQ simulation is found to be smoother at low 
redshift (z < 3) than at high redshift (z > 3). The final mass of 
the dark matter halo in the NutCO run is M v i r ~ 4 x 10 11 Mq 
(see Appendix for a time evolution of the halo mass), which is 
thought to be close to the mass at whic h gas accretion transitions 
from the cold mo de to the hot mode JBirnboim & Deke] 120031 : 
lOcvirk et alj|2008h . Since the NutCO run does not include super- 
nova feedback, satellite galaxies can retain more gas than in the 
NutFB run, and thus the resulting relative contribution from the 
cold, dense filamentary gas to the total angular momentum gets 
slightly smaller at high redshift than the feedback run (Fig. El- 
However, there is little difference in / gas ,out or /dm between these 
two runs for z > 3. Fig. El also shows that at z ^ 2, the diffuse 
hot gas phase becomes the dominant reservoir of angular momen- 
tum in the region (0.1 < r/7? v ir < 1). Considering that the post- 
shock temperature for an isothermal sphere of mass a few times 
10 11 Mq (our halo reaches 2 x 10 11 M W at z ~ 2, c.f . Appe ndix A) 
is Tshocked ~ 3T v i r /8 ~ 10* K JPekel & Birnboimll200(ih . where 
T v i r (— 35. 9x (I4/[km • s -1 ]) 2 ) is the virial temperature, and that 
the run does not feature any SN feedback, we can safely identify the 
halo hot gas at z ^ 2 with material which has been shock-heated. 


1 Run with identical initial conditions and mass resolution to the NutFB 
run discussed in the previous subsection, but with lower (48 pc) spatial res- 
olution and no supernova feedback 
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Figure 5. Same as Fig. [3 but for the NutCO run. Different colours denote different components: dark matter (thick grey), gas at 0.1 < r / R./ u < 1 (thick 
blue), central gas (r/f? v ; r < 0.1, thin blue), and stars (orange). Note that the dark matter halo and the halo gas includes the contribution from satellites, 
while they are excluded for the measurement of J of stars. We also display j of all baryons (star+gas) as a red dot-dashed line. Specific angular momentum 
corresponding to A' = 0.04 is included as a grey dashed line. At low redshift, the misalignment between the central gas and each component is small, and gas 
in the outer region shows larger angular momentum than dark matter in both panels. 


Nonetheless, as Fig.[5]cleariy demonstrates, this hot halo gas within 
0.1 < r / .Rvir < 1 still has a specific angular momentum higher 
than that of the dark matter, indicating that shock-heating of the ac- 
creted cold gas, even when it takes place as early as z = 2, does not 
erase the difference in specific angular momentum of dark matter 
and gas. 

As is the case in the NutFB run, we find that jbar is compa- 
rable to jdm up to z=0, whilst y gas, out is always greater than jd m . 
Once again, given that gas accretion at the virial radius of the halo 
is never perfectly coplanar with the galactic disc, we expect the 
(small) off-axis components (J x and j y ) to vanish before gas be- 
comes part of the central galaxy, and thus the modulus of gas spe- 
cific angular momentum not to be conserved in general. We also 
expect gas accretion history to reduce the modulus and smooth the 
time evolution of j ga s,cenwith respect to y ga s,out. Fig.[7]shows the 
radial profiles of the j distributions of gas and dark matter within 
the virial radius of the dark matter halo, normalised to their respec- 
tive values of j at 0.1R v i r . These profiles are split in a low and a 
high redshift bin (z < 3 and 2 > 3) and stacked. Two striking fea- 
tures emerge from the analysis of Fig. [7] First the dependence of 
the shape of all the specific angular momentum profiles on redshift 
is minor. Second, while the j of the dark matter component drops 
almost proportionally with radius, the j profile of the gas is flat in 
the outer parts of the halo (r > 0.1 R v i r ). Although we did not 
included the curve in this paper, this is in stark contrast with the 
behaviour of the gas specific momentum in the NutAD run which 
closely follow the j profile of the dark matter halo. This implies 
that as long as the gas can radiatively cool, its specific angular mo- 
mentum is conserved whilst it is being advected to the inner region 
(r < 0.1 .Rvir), regardless of whether shock heating occurs or not. 
Therefore, any redistribution of angular momentum must happen 
within the central region of the halo, close to the galactic disc. The 
detailed analysis of this complex process is beyond the scope of 
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Figure 6. Same as Fig. [4] but for the NutCO run. The hot gas phase (orange 
lines) begins to develop at z ~ 3, and accounts for almost all of the angular 
momentum in 0.1 < r/R v - 1T < 1 by z = 0. 


the current paper and will be presented elsewhere (Tillson et al. in 
prep.). 
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Figure 7. Stacked distributions of normalised specific angular momentum 
as a function of radius from the NutCO run. The specific angular momentum 
of gas and dark matter is shown as solid and dashed lines. Blue lines display 
j at high redshift (z > .3), while orange lines show j at low redshift (z < 3). 
Each point corresponds to the median, and interquartile ranges are included 
as error bars. The normalised j of gas is flat at 0.1 1$ r/R v ; r T 1, and gets 
diminished in the region where the galactic disc is located, suggesting that 
j of gas is reasonably conserved up to r = 0.1 R v \ r . 

3.4 Generalisation of the results 

Now that we have established that gas in a Milky Way-like halo has 
a larger specific angular momentum than dark matter, it is vital to 
know whether this property is generic or not. To build a statistically 
representative sample of haloes, we make use of three large cosmo- 
logical hydrodynamics simulations called HoRizoN-MareNostrum, 
Cosmo25, and Cosmo50. The HoRizoN-MareNostrum simulation 
contains 326, 1921, 3307, and 4015 haloes with M v i r > 10 11 Mg 
at z=6.1, 3.8, 2.5, and 1.5. Since the simulation has only been car- 
ried out down to z= 1.5, we use the Cosmo25 and Cosmo50 run to 
supplement the sample at lower redshifts. The number of haloes 
with M v ir > 1O 11 M 0 in the Cosm25 run is 378 and 459 at z=1.5 
and z=0 respectively. For the Cosmo50 run, we restrict our analysis 
to 481 more massive haloes (M v i r > 10 12 Mq) at z=0 because of 
the poorer resolution of the simulation. 

The middle panel of Fig. [8] shows the median and the in- 
terquartile range of the spin parameter distribution of the dark mat- 
ter haloes. As is well known in the literature, it shows that the typ- 
ical value for the spin is A ~ 0.04 with very little dependence 
on halo mass and redshift. On the other hand, the spin parameter 
of gas in the HoRizoN-MareNostrum simulation (filled circles) is 
2-4 times higher in general than that of the dark matter halo and 
clearly depends on halo mass. The ratio turns out to be even more 
significant at z=0 in the Cosmo25 (asterisks) and Cosmo50 run 
(diamond), demonstrating that our finding that gas has larger spe- 
cific angular momentum extends to all halos over the entire redshift 
range. Moreover, it also substantiates our claim that gas brought in 
later in the life of the halo has more specific angular momentum be- 
cause of the large scale origin of this angular momentum which has 
been investigated in a companion paper IPichon et al.ll201 il l. Note 



Figure 8. Spin parameter distributions of the gas within 0.1 < r/R v i r < 1 
(top), the dark matter halo (middle), and the ratio of the two (bottom). Filled 
circles correspond to the results from the HORIZON-MareNostrum simula- 
tion, and asterisks and diamonds show the results from the Cosmo25 and 
Cosmo50 run, respectively. Different colours denote different redshifts as 
indicated on the figure. The results from the NutCO run are also included 
as stars with the same colour-coding. Solid lines show linear fits to the sim- 
ulated data at different redshifts as given by Equation 0 - Statistically, the 
specific angular momentum of the gas within 0.1 < r / R vlr < 1 is at least 
twice that of the dark matter halo host, at all redshifts. 

that since the HoRizoN-MareNostrum and Cosmo25 simulations 
are run using different cosmologies (WMAP1 and WMAP5 respec- 
tively) and feature different mass resolutions, we included data at a 
common redshift (z = 1.5) to ensure that the differences between 
the z=1.5 HoRizoN-MareNostrum and the z=0 Cosmo25 halos does 
not arise from a different choice of cosmology/simulation parame- 
ters. 

For practical purposes, we compute a redshift-dependent fit to 
the specific angular momentum gas to dark matter ratio for halos 
with M v i r < 10 13 Mq 
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Figure 9. True colour composite image of the NutCO galaxy in the U, B 
and V at z = 0. The image measures 15 kpc on a side. Note that internal 
extinction by dust is not accounted for in this image. Notice the presence of 
a large bulge (I-band B/D ratio 0.78) typical of a galaxy with morphological 
type Sa/b (see text for detail). 


ga , s,OUt ^ max[3.5+24.3(l+z)-°- 8 -1.97(l+z)“°- 7 logM vir , 2].(3) 

^DM 

We emphasize that all the extra specific angular momentum 
brought in by the gas will not be devoted to spinning up the central 
galactic disc, given that its typical misalignment with the central 
gaseous disc is measured to be « 40° for intermediate mass ha- 
los (M v ir ~ 1O 11 M0) in the HoRizoN-MareNostrum simulation. 
Moreover, in the vicinity of the central disc, hydrodynamical in- 
teractions with the circum-galactic gas will become important and 
may substantially redistribute angular momentum. A careful nu- 
merical investigation of angular momentum advection in the cen- 
tral region is therefore needed to determine whether or not a simple 
semi-analytic approach is capable of correctly describing the evo- 
lution of angular momentum of disc galaxies using the new initial 
conditions we provide with this fit. 

We also note that limited spatial resolution of our large volume 
cosmological simulations will lead to an artificial increase of the 
gas angular momentum, especially in low mass halos. Indeed, the 
softening of the gravitational force produces very extended discs 
at the centre of these halos, which would probably be contained 
within 0.1 R vlr at higher resolution. However, comparing the spin 
of the gas in the NutCO halo (large empty stars in Fig.[8]t at various 
redshifts with that measured for a sample of halos of comparable 
mass available in these cosmological simulations (solid circles and 
asterisks in Fig. [8}, we conclude that resolution effects most likely 
account for a minor fraction of the gas spin (~ 20%). 


3.5 Comparison with observations 

Ultimately, to decide whether or not our simulations yield realistic 
galaxies, we have to compare their properties with observations. 
Unfortunately, with the high resolution we need to properly address 
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Figure 10. Comparison between the NutCO galaxy a nd th e observational 
bright disc galaxy sample described in lKassin et alJ UOQbl) . Observational 
data is in grey (lower and upper estimates of j are given by empty and solid 
circles respectively, see Kassin et al. 201 1 for details). For the simulation 
data, different symbols (circles and right arrows) stand for different esti- 
mates of the velocity of the galaxy (circular velocity and maximum rotation 
velocity). We also include estimates of V 2 . 2 , the velocity measured at 2.2 
times the exponential scale length in I band, for z < 1. Note that the simu- 
lated galaxy displays a dense large bulge which overgrows, at least partially, 
because of the lack of modelling of any kind of internal feedback (radiation, 
stellar winds, supernovae, cosmic rays, etc.) in the NutCO run. 



the issues related to angular momentum, cooling and feedback, we 
can only perform such a comparison at z = 0 with one of our 
runs (NutCO) at the time being. As this simulation does not model 
any feedback mechanism, by z=0 the simulated galaxy has formed 
too many stars (M star = 6.8 x 10 10 Mq). At z = 0, the /-band 
bulge-to-disc ratio (B / D) of the central galaxy in the NutCO run 
is ~ 0.78 (see Fig.[9j. Such a value is close to the typical B/D of 
Sa/Sb type spirals, but the disc scal e lengt h (2.0 kpc) turns out to 
be smaller than those observed l iGraham & Worlevll2008l) . 

Nevertheless, bearing this caveat in mind, we present in 
Fig. 1101 measurements of the j versus V relation for the NutCO 
galaxy at various redshifts, com pared to the observ ational bright 
disc galaxy sample gathered by iKassin et al] J2006I) (at z ~ 0). 
Kassin et al. derived the specific angular momenta for these galax- 
ies from Halpha + HI r otation cur ves and radial stellar mass distri- 
butions given in ( Kassin et al1 l2006l ). The average errors in j mea - 
sured from the data are « 60kpckms _1 . As is done in the ob- 
servations, both the amount of specific angular momentum and the 
velocity of the simulated galaxy (Umax, maximum of the rotation 
curve) are estimated using the stellar component, except for the 
filled circle symbols where the velocity is measured as the circular 
velocity at the virial radius of the dark matter halo. For a fair com- 
parison with observations, we also include velocities measured at 
2.2 times the disc scale length in the / band (V 2 . 2 ) at z < 1. 

Looking at Fig. [TO] one clearly sees that the velocity estimate 
plays a crucial role in our ability to assess whether simulated galaxy 
stars have the correct amount of angular momentum. This is be- 
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Figure 11. Accretion- weighted specific angular momentum for dark matter (black solid line) and gas (blue solid line) at the virial radius in the NutFB (left 
panel) and the NutCO (right panel) runs. The specific angular momenta of the dark matter and the halo gas (jgas,out) shown in Fig.[3]and[5]are included as 
grey dashed and blue dashed lines. Note the similar evolution of the specific angular momentum of the halo gas and of both accreted gas and dark matter. 


cause the observed slope of the relation (« 2) is quite large so any 
error on the velocity will translate into a much larger error on the 
specific angular momentum. More specifically, our simulation suf- 
fers from forming too many stars and as a result any velocity esti- 
mated in the central region of the halo is bound to be too large. If, 
on the other hand, we use the circular velocity at the virial radius 
of the h alo to b y pass the pro blem (note that this is current practice, 
see e.g. iDutton & van den Boscbll2009l ). we find that the level of 
angular momentum of the simulated stars is in fair agreement with 
the observations (Fig. ED. We point out that this is a consequence 
of the lossless transport of a large amount of specific angular mo- 
mentum by gas from super halo scales right into the inner regions, 
followed by its redistribution in the vicinity of the disc, as discussed 
in the previous section. As a result of this effect and the low effi- 
ciency of the star formation process on disc scales, stars end up 
with a factor 2 to 3 less angular momentum than their dark mat- 
ter halo counterpart which is very close to the discrepancy between 
dark matter only simulation s and observations, as pointed out by 
INavarro & Steinmeta 0000 ) and Kassin et al. (201 1 , in prep.). 

Finally, it is interesting to note that when we use the V max es- 
timates for the velocity of the galaxy, although we do not match 
the zero point of the relation, simulated data points at different red- 
shifts move along the observed sequence, suggesting that we would 
predict very little time evolution of the sequence. On the contrary, 
V c estimates suggest there exists a strong evolution of the relation 
with redshift as the flat velocity profile of the dark matter halo is 
already in place very early on. 


4 ORIGIN OF THE DISCREPANCY BETWEEN DARK 
MATTER AND GAS SPECIFIC ANGULAR MOMENTA 

Having robustly demonstrated that the amount of specific angular 
momentum carried by gas and dark matter differ throughout the 
virialised halo and at all redshifts, we now turn to identifying the 


cause of such a difference in j. There are essentially three possi- 
bilities: i) the gas carries a larger amount of specific angular mo- 
mentum than the dark matter at the time of accretion, or ii) the dark 
matter component does not conserve j as well as the gas or iii) the 
difference in j reflects a different spatial re-distribution of the an- 
gular momentum of gas and dark matter within the virialised halo. 

In order to test the first possibility, i.e. whether the gas already 
carries a larger amount of specific angular momentum than the dark 
matter when it is accreted by the halo, we compute the accretion- 
weighted specific angular momentum modulus of both components 
at the virial radius as 

_ X], miwO(-'iv ) Vi X Vi 

3 Vr miVrO(-Vr) 

where 0 is the Heaviside step function. This is done by calculating 
the contribution of infalling ( v r < 0) gas or DM particles located 
within the thin shell defined by 0.95 < r/R v i r < 1.05. Since 
the collisionless dark matter particles can potentially get accreted 
several times as they come in and out of the halo, we track each 
particle individually and flag them the first time they cross the virial 
sphere so as not to re-accrete them later. For the collisional gas, this 
problem does not arise as it remains in the center of the halo unless 
blown out/evaporated by feedback or shock-heating, so a simple 
radial velocity cut suffices. 

Fig- ED shows that the specific angular momenta of the ac- 
creted gas and dark matter agree reasonably well both in the NutFB 
and NutCO run, at all times. Although we have not included it in 
this paper, this turns out to be the same in the NutAD ran. It is 
noteworthy that both freshly accreted gas and dark matter bring 
in a larger amount of specific angular momentum than that of the 
dark matter halo, j'd m , which is the common assumption used in 
(semi-)analytic studies. The fact that j gas , out carries nearly the 
same large amount of specific angular momentum as (j) sup- 
ports the view that the angular momentum of the gas accreted by 
the central galaxy is acquired before it enters the dark matter halo. 
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As argued elsewhere JPichon et all 1201 II) . this favours the large- 
scale motion of the cosmic web as the most likely origin of the 
angular momentum of galaxies. 

The second possibility, i.e. that the amplitude of the dark mat- 
ter component specific angular momentum is conserved to a lesser 
degree than the gas, is easily ruled out on the grounds that, as ar- 
gued in Section [3721 when the time evolution of the total specific 
angular momentum of all the baryons is calculated, it matches that 
°f jdm quite accurately. It seems therefore impossible to argue that 
dark matter loses more specific angular momentum than the gas as 
the main reason why j gas , out is larger than jdm- 

We are then left with the third possibility: the dissipative na- 
ture of the gas leads to a transfer of specific angular momentum dis- 
tinct from that of the dark matter. As previously mentioned, cold, 
dense filamentary gas flows into the inner part of the dark matter 
halo in a dynamical time without und ergoing shock-heating (e.g. 
iBrooks et alj|2009l : IPowell et al .11201 ll) . The gas subsequently set- 
tles onto the central gaseous disc, and, as a result, the outer region 
(0.1 < r / Rvii < 1) is occupied by gas recently accreted, which 
has a large amount of specific angular momentum. On the other 
hand, dark matter particles pass through the central region of the 
halo and depending on their exact orbital properties, can end up 
populating the outer regions despite having entered the halo at ear- 
lier times with smaller angular momentum. This is shown in Fig. ED 
where we plot the angular momentum as a function of position for 
all the DM particles belonging to the NutCO halo at z=0, colour 
encoding the time at which they where accreted. Whilst a lot of par- 
ticles in the outer regions were accreted fairly recently (less than 5 
Gyr ago), there exists a non- negligible fraction of them which were 
accreted at much earlier epochs (more than 7 Gyrs ago). Thus, de- 
spite the fact that dark matter first falls into the halo with a similar 
amount of specific angular momentum as the gas (Fig. ED, mixing 
with ‘older’ particles makes jdm become systematically smaller 
than y gas ,out. Note that this also explains the similar evolution of 
the specific angular momentum of the halo gas Ogas.out) and the 
newly accreted gas at high redshift (Fig. ED- 

Obviously the cold, dense filam entary gas accre- 

tion disappears at low redshift (e.g. jKimm et al J 1201 ll : 
iFaucher-Giguere & Kereslholllfste w art et all20f(ih . and therefore 
one expects that the high redshift segregation of specific angular 
momentum between gas and dark matter previously advocated 
becomes weaker as time elapses. Indeed, in the extreme case 
where gas is prevented from cooling (NutAD run: Fig. [2} the 
evolution of its specific angular momentum closely follows that of 
the dark matter, indicative of a similar degree of mixing of low and 
high specific angular momentum material for both components. 
However, when gas can radiatively cool, one expects this mixing 
process to be much less efficient as it will generally take place 
on timescales (typically a few sound crossing times of the halo) 
which are longer than the cooling time of the hot central gas. 
This means that the shock-heated material at the centre of the 
halo, which has a lower specific angular momentum than that in 
the outer regions, is cooled and accreted onto the central disk 
before it can mix with higher specific angular momentum gas. 
This creates a ‘cooling flow’ towards the central galaxy and as a 
result preserves the distribution of gas specific angular momentum 
that was set in at high redshift (Figs. [5] and |7J. Note that this 
also explains the increasing discrepancy between dark matter and 
gas spin as the mass of the halo drops that we measured in large 
scale cosmological simulations (Section 13.41 Fig. GQ: as the halos 
become larger, cooling timescales increase (because of higher gas 
temperatures and smaller central densities) allowing more and 



0.0 0.2 0.4 0.6 0.8 1.0 

r/r vir (z=0) 


Figure 12. Distribution of angular momentum of DM particles as a function 
of their position for the virialised DM halo in the NutCO run at z = 0. Col- 
ors indicate the lookback time at which particles were accreted for the first 
time. Dark and light grey contours enclose 50 and 90 percent of the mass 
respectively. Note the large range of angular momentum values spanned by 
individual particles at any radius. 

more mixing between low and high specific angular momentum 
gas to occur. As a result, the gap between j ga s,out and 7dm is 
reduced for the most massive objects. 

The redshift dependence of the spin parameter of the gas can 
also be understood in this framework. Again, since the cooling per- 
mits the gas to flow into the centre, 7g as ,out tends to reflect j of 
newly accreted gas. This means that, for any dark matter halo of 
a given mass, its gas content at lower redshift will be composed 
of material carrying larger specific angular momentum. In other 
words, the redshift dependence of gas spin shown in Fig. [8] is un- 
derstood as a large scale structure driven increase in the specific 
angular momentum of freshly accreted gas. 

Even though we do not investigate what causes this behaviour 
in this paper, we remark that the systematic discrepancy between 
jdm and (j) measured in the NutCO run at all times means that 
the dark matter component of the halo (and its baryons as a whole) 
loses somewhere between half and two thirds of its specific angu- 
lar momentum amplitude, depending on redshift (Fig. ED . We note 
that even thou gh thi s result seems in agreement with the measures 
of lBook et al .1 120111) for a handful of well resolved dark matter ha- 
los in pure N-body simulations, we do not find, like these authors, 
that it can be attribute d to e xternal torques, following the original 
argument of lPeebleD l l 19691) . The reason for this disagreement lies 
in the continuous nature of the specific angular momentum loss of 
the dark matter component. As Fig.QT]demonstrates, the accretion 
weighted specific angular momentum of dark matter is higher than 
the specific angular momentum of the dark matter halo at all red- 
shifts, regardless of whether we are in a phase where tidal torque 
theory predicts the angular momentum of the halo should grow or 
not. We are therefore pushed to conclude that the angular momen- 
tum Toss’ suffered by the halo has to be attributed primarily to a 
‘vector cancellation’ effect: particles with similar amounts of angu- 
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lar momentum (modulus) but pointing away in opposite directions 
contribute only a small amount to th e ‘net’ a ngular momentum of 
the virialised halo (e.g. Ivitvitska et al.ll2002h . Note that the accre- 
tion history effect will also play a role in reducing the modulus of 
jdm but cannot account for all the decrease. 

Finally, the admittedly low spatial resolution (1-2 kpc) of the 
large scale cosmological hydrodynamics simulations we use to gen- 
eralise our results, prevents dark matter (and gas) from forming 
more concentrated structures. As a result, galaxies in small halos 
(JVdm ~ 2 x 10 3 ) show a more extended disc than their Nut se- 
ries counterparts. However, we believe that the properties of halos 
we discussed in this paper, since they are integrated over a large 
fraction of the halo radius, are robust. A caveat to bear in mind 
however, is that the physical processes that can prevent the cool- 
ing catastrophe in cluster environments are likely to change the 
ratio Agas.out/Adm- For example, as extensively discussed in the 
literature, feedback from active galactic nuclei is thought to stir/ 
heat a large amount of cooling gas, redistributing it in the outer 
parts of the halo. As previously argued, such a mixing mechanism 
will contribute to narrow the gap between gas and dark matter spin. 
Such a potent feedback mechanism is not included in the cosmo- 
logical simulations presented here, and thus a larger fraction of 
baryons than observed are converted into stars in massive haloes 
(Dubois e t alj201(il) . This should be taken into account before con- 
cluding that the ratio (A gas ,out/Adm) for massive haloes has con- 
verged to ~ 2 as Fig. [8] (bottom panel) advocates. 


5 CONCLUSIONS 

In this paper, using high-resolution cosmological hydrodynamics 
re-simulations of individual galaxies and large scale cosmological 
hydrodynamics simulations, we have revisited the standard theory 
for angular momentum evolution of gas and dark matter within viri- 
alised structures. 

According to this standard picture, gas acquires angular mo- 
mentum through nearby tidal fields in the same way as dark mat- 
ter and is shock-heated to the virial temperature during the col- 
lapse of the halo. This gas then radiatively cools and settles into 
a rotationally supported disc, conserving its specific angular mo- 
mentum in the proce ss dFall & E fstathiou 1 98(1 ; iDalcanton et al.l 
1 19971 : InIo et al.lll998h . Such assumptions naturally lead to predict 
that dark matter and halo gas share the same specific angular mo- 
mentum. Investigating the issue, our main conclusions can be sum- 
marised as follows: 

• Indeed, the time evolution of the halo gas specific angular mo- 
mentum is quite comparable to that of the dark matter, but only 
when gas is not permitted to radiatively cool, i.e. it is artificially 
forced to shock- heat during the collapse (NutAD run). 

• When radiative cooling is turned on, the angular momentum 
transport of gas becomes distinctively different from that of dark 
matter (NutCO run). Even though gas and dark matter still bring 
in a similar specific angular momentum ({j) v ) when they are ac- 
creted, this amount turns out to be systematically larger than the 
specific angular momentum of the virialised dark matter halo as a 
whole (jdm) b y a factor 2 to 6, depending on redshift and halo mass 
(see also lStewart et alj201 ll) . 

• Regardless of whether cooling is turned on or not, the modulus 
of the specific angular momentum inside the virial sphere is not 
conserved for dark matter or gas. The amount lost ranges between 
one half and two thirds depending on redshift. We attribute this 
continuous Toss’ of angular momentum to a ‘vector cancellation’ 


effect (the angular momentum of newly accreted material is never 
perfectly aligned with the angular momentum of the whole halo) 
augmented by an accretion history effect (material accreted earlier 
carries less angular momentum so it weighs the average modulus 
of the specific angular momentum down at any epoch). 

• When radiative cooling is turned on, provided one compares 
the dark matter with the total baryon (gas and stars) specific angu- 
lar momentum, their amplitudes are similar at all times. In other 
words, we find no evidence of baryons losing more (or less) an- 
gular momentum than dark matter within the virialised halo at any 
epoch. 

More specifically, our analysis reveals that at high redshift, 
the discrepancy between gas and dark matter specific angular mo- 
mentum within 0.1 < r/R v ir < 1.0 (j gas ,out and jd m , respec- 
tively) arises because the angular momentum-rich freshly accreted 
gas flows into the central region through cold, dense filamentary ac- 
cretion without being redistributed throughout the halo by shock- 
heating or supernovae feedback (NutFB run). As the dark matter 
halo grows at lower redshifts, a progressively larger amount of gas 
undergoes shock-heating, but since the central dense gas can still 
cool and collapse onto the disc on timescales shorter than a halo 
sound crossing time, the difference between j gas , out and jdm is, 
by and large, preserved. This competition between central cooling 
and redistribution of gas through shock heating and/or feedback in- 
duces a mass dependence of the discrepancy between gas and dark 
matter specific angular momentum: it is reduced for more massive 
halos where shock heating dominates. However, for a halo of fixed 
virial mass, this discrepancy is also larger at low redshifts than it is 
at high redshifts because the amount of specific momentum carried 
by the newly accreted material grows faster than the halo average 
with time (cf. Pi chon et ah . J201 11) for the origin of this effect). This 
generic behaviour is encapsulated in Equation 0. which provides 
a fit to our simulated data. 

Our efforts to probe the spatial distribution of angular momen- 
tum within the virialised halo also led us to define an inner region 
(r/R v i r < 0.1). Gas in this region, j gas , C en, generally has more 
specific angular momentum (by a factor 2 or so) than the dark mat- 
ter halo as a whole, and only stars have less (also by about a factor 
2). We attribute this dichotomy within the baryon component to the 
long global star formation timescale of the disc compared to the 
gas accretion time scale, driving stars to preferentially form from 
gas accreted at a mu ch earlier stag e wit h low er angular momen- 
tum (see also iDutton & van den Boscbl j2009i) who reach similar 
conclusions with a SAM). As a result, a comparison with the data 
at z=0 reveals that the level of specific angular momentum of our 
high resolution simulated stellar disc matches quite well that of the 
observations. However the difficulty of obtaining realistic galaxy 
rotation curves still prevents us from satisfactorily reproducing the 
observed j versus V relation (see Fig. ED. Part of this failure is to 
blame on the fact that the only re-simulation we have been able to 
run down to z=0 so far (NutCO run) does not include any modelling 
of feedbac k (stellar winds, supemovae, cosmic ray, etc . ) (as pointed 
out b y e.g. iGovemato et alj|2007l : lAgertz et~aT1l20 1 ll : iBrook et all 
1201 ll) . For this reason, too many baryons (the universal fraction 
Hs/!l ra ) remain in the host halo and concentrate in the galaxy, po- 
tentially affecting the morphology of our galaxy which features a 
large bulge (morphological type Sa/Sb) but in any case leading us 
to overestimate disc rotational velocities. 

Overall, our results demonstrate the need for the standard pic- 
ture of galaxy disc formation and evolution to undergo a major 
overhaul as the fundamental hypothesis upon which it is based 
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(specific angular momentum of the gas and dark matter within 
the virialised halo are the same, gas specific angular momentum 
is conserved) seem to be incorrect. This should have important 
consequences for galaxy disc sizes and therefore star formation 
timescales estimates used in virtually all galaxy formation mod- 
els, although it is possible (albeit very unlikely) that errors induced 
by these two assumptions exactly cancel out. The new picture we 
advocate should take into account how gaseous cold flows effec- 
tively carry larger than previously thought amounts of angular mo- 
mentum originating from the large-scale motion of the cosmic web 
down to the central region of virialised structures. The fundamen- 
tal role played by gas cooling at low redshift which, for galaxy size 
halos, preserves this transport of angular momentum in spite of the 
presence of a pervasive shock-heated corona which entirely fills the 
virialised halo, should be described. Finally, since gas accretion at 
the virial radius is not perfectly aligned with the axis of rotation 
of the galaxy, some angular momentum will inevitably cancel out 
in the vicinity of the disc. Careful numerical investigation of how 
angular momentum is transported and/or cancelled in the central re- 
gion of the halo should help shed light into this fundamental aspect 
of disc formation theory. 
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APPENDIX: MASS ACCRETION HISTORIES 

In the main body of the text, we show the time evolution of specific angu- 
lar momentum for the various components of a Milky Way-like DM halo 
(Figs. [5] an d[5). In particular, in Section 3.2, we argue that the difference in 
angular momentum between the central gas and the gas in the outer region 
of the halo arises because the mass of the central gas, which is composed of 
low angular momentum gas accreted at earlier time, is not negligible com- 
pared to the mass of accreting gas, and therefore reflects the history of gas 
accretion. To substantiate this claim, we include the mass assembly history 
of the various components in Fig. m As can be seen, the amount of gas in 
the central region (M gas — Mg as (r > 0.1_R v ir)) is comparable to the gas 
mass in the outer region (r > 0.1i? v i r ) at all time. Thus, even if the newly 
accreted gas carries a larger amount of specific angular momentum, one 
expects that there will be a non-negligible time-delay for the entire central 
region to be spun at the same level. Finally, note that by multiplying Fig. ca 
with Figs.|2and[5] one can recover the total amount of angular momentum 
of each component. 




Figure 13. Mass assembly histories of the various components in the NutFB 
run (top) and the NutCO run (bottom). Stars, dark matter, and gas in satellite 
galaxies are included in the measurement. Note that the amount of gas mass 
in the central region (M ga s — M gas (r > 0.1_R v ir)) is comparable (a factor 
2-3 lower) to the gas mass in the outer region at all times. 
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